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Abstract 

Colorectal cancers are the third most common type of cancer. They 
originate from intestinal crypts, glands that descend from the intestinal 
lumen into the underlying connective tissue. Normal crypts are thought 
to exist in a dynamic equilibrium where the rate of cell production at the 
base of a crypt is matched by that of loss at the top. Understanding how 
genetic alterations accumulate and proceed to disrupt this dynamic equi¬ 
librium is fundamental to understanding the origins of colorectal cancer. 
Colorectal cancer emerges from the interaction of biological processes that 
span several spatial scales, from mutations that cause inappropriate intra¬ 
cellular responses to changes at the cell/tissue level, such as uncontrolled 
proliferation and altered motility and adhesion. Multiscale mathematical 
modelling can provide insight into the spatiotemporal organisation of such 
a complex, highly regulated and dynamic system. Moreover, the afore¬ 
mentioned challenges are inherent to the multiscale modelling of biological 
tissue more generally. In this review we describe the mathematical ap¬ 
proaches that have been applied to investigate multiscale aspects of crypt 
behaviour, highlighting a number of model predictions that have since 
been validated experimentally. We also discuss some of the key mathe¬ 
matical and computational challenges associated with the multiscale mod¬ 
elling approach. We conclude by discussing recent efforts to derive coarse¬ 
grained descriptions of such models, which may offer one way of reducing 
the computational cost of simulation by leveraging well-established tools 
of mathematical analysis to address key problems in multiscale modelling. 

Keywords: OfF-lattice models; continuum models; multiscale models; coarse- 
graining; intestine; colorectal cancer. 
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1 Introduction 


The past decade has witnessed remarkable progress in experimental studies of 
intestinal tissue self-renewal and colorectal carcinogenesis. This is due in part 
to increased molecular understanding of intestinal cell biology, as well as ad¬ 
vances in techniques such as clonal analysis, [smsiiM] in vitro tissue culture [55] 
and intravital microscopy. [53] Mathematical modelling offers a complementary 
tool with which to unravel the complex interactions between processes at the 
intracellular, cellular and tissue scales that underlie the spread and fixation of 
mutations during the first stages of colorectal cancer. Models can be used to 
develop abstract representations of biological systems, test competing hypothe¬ 
ses, and generate new predictions that can then be validated experimentally. 
To this end, a variety of different individual cell-based and multiscale modelling 
approaches have recently been developed for studying how processes at the 
level of a single cell affect the tissue-level behaviour of the intestinal epithelium. 
This review aims to highlight the mathematical challenges inherent in multi¬ 
scale modelling of intestinal self-renewal and colorectal cancer, which serves as 
a representative example of the more general study of stem cell dynamics, tissue 
homeostasis and carcinogenesis. 

The remainder of the paper is structured as follows. In Section |2| we sum¬ 
marise the relevant background biology. In Section |3| we briefly discuss the 
variety of mathematical modelling approaches that have been applied to the 
study of colorectal cancer, placing multiscale models into a wider context. In 
Section 0] we review multiscale modelling efforts, focusing on the competing 
approaches and on the model predictions that have since been validated experi¬ 
mentally. In Section [5] we highlight some of the key mathematical and computa¬ 
tional challenges associated with the multiscale modelling approach and discuss 
recent efforts to derive coarse-grained descriptions of such models, which aim 
to improve analytic tractability while retaining the essential features of the un¬ 
derlying models. Finally, in Section |6| we offer concluding remarks on future 
directions for this work. 


2 Background biology 

Colorectal cancers initiate in the epithelial tissue that covers the luminal sur¬ 
face of the intestinal tract. The rapid renewal of this epithelium involves a co¬ 
ordinated programme of cell proliferation, migration and differentiation, which 
begins in the tiny intestinal crypts, finger-like invaginations of the intestinal 
mucosa lined by a monolayer of epithelial cells (Fig. ITal) . Crypts function to 
absorb water and exchange electrolytes from the faeces, and to produce mucus 
to lubricate faeces as they move through the large intestine (colon). Colorectal 
carcinogenesis occurs as a consequence of changes that disrupt normal crypt 
dynamics. Thus, the tissue architecture and dynamics of intestinal crypts are 
central to the origins of colorectal cancer. |M) 

A population of stem cells located at the base of each crypt divide frequently 
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Figure 1: (a) Scanning electron micrograph of murine colonic epithelium, show¬ 
ing the structural crypt organization, (b) Schematic of the spatial organization 
of the colonic crypt. Stems cells (red) residing at the crypt base generate rapidly 
proliferating transit amplifying cells (green), which ascend the crypt and differ¬ 
entiate into the mature cell lineages (white) occupying the surface epithelium 
(in the colon) or adjacent villi (in the small intestine), (c) Schematic of the 
process of monoclonal conversion, where the progeny of a single stem cell (blue) 
replaces all other cells in the crypt as a result of neutral drift, (d) Schematic of 
the process of crypt Hssion, where excess proliferation at the crypt base results 
in buckling, basal bifurcation and longitudinal division; this process facilitates 
the spread of mutations beyond individual crypts. Image (a) reproduced with 
permission from Magney et al.[46] ©Wiley (1986). 
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to generate transit amplifying cells that migrate up the crypt wall, [5] as shown 
in Fig. [Tbl As they ascend, these cells undergo a sequence of rapid divisions 
before they stop dividing and differentiate terminally into specialised progeny, 
such as colonocytes, enteroendocrine cells or Goblet cells, that occupy the upper 
portion of the crypt. mm In the colon, cells at the crypt orifice either undergo 
apoptosis or are shed into the lumen and transported away, while in the small 
intestine they move onto an adjacent villus structure. Crypts may be generated 
by a process called crypt fission, in which a crypt bifurcates basally and then 
divides longitudinally. Crypt fission occurs sporadically in the adult intestine 
as well as in response to epithelial damage. m Crypt extinction is also possible 
and occurs if all actual and potential stem cells are eliminated, for example after 
exposure to radiation or cytotoxic drugs. m 

Intestinal stem cells are in continual neutral competition with one another to 
retain a privileged position within the niche at the crypt base. [S] Each stem cell 
division results, on average, in the loss and replacement of an individual stem 
cell lineage. 0 IMl Sil IBH] A consequence of this competition is that over time, 
the progeny of a single stem cell within a crypt can take over the entire crypt. 
This process is termed monoclonal conversion, since the resulting crypt consists 
of a single clonal population. Mutant monoclonal crypts are likely to constitute 
the earliest stage of colorectal adenomas, which then expand further through 
crypt fission rather than unconstrained growth of cells per se. [SnilSlElEH] The 
rate of fission is significantly altered in such mutant crypts; for example, the 
basal rate of crypt fission in the mouse colon is increase 30-fold by the oncogenic 
KRAS mutation. [69] This is likely due to the increased stress on crypt walls due 
to the persistence and excess proliferation of mutant cells, which can cause crypt 
buckling and fission. 

A number of biochemical signalling pathways have been implicated in the 
regulation of intestinal stem cell and crypt dynamics. Of particular importance 
in the context of colorectal cancer is Wnt signalling, which is initiated when 
extracellular Wnt factors bind to specific receptors on the cell surface. This 
triggers a cascade of intracellular reactions that leads, via the transcription of 
Wnt target genes, to the regulation of proteins involved in cell-cycle control, 
migration and apoptosis. Mutations in key components of the Wnt pathway, 
such as APC and /3-catenin, have been shown to be the first step in colorectal 
carcinogenesis in the majority of colorectal cancers, including most hereditary 
cases. [35] Notably, the stem cell niche in the colonic crypt is regulated by a 
spatial gradient of extracellular Wnt factors along the crypt axis. [^ [74] 

Despite the central role of stem cells in maintaining the integrity of the 
intestinal epithelium, for many years their precise number and location within 
the crypt has remained unclear, due to a lack of unique marker genes and the 
absence of stem-cell assays. As the transgenic approach used successfully in 
mice [aSTlITT] is impractical for human studies, it has remained uncertain if 
and how the biology of the human intestinal crypt is mirrored by its murine 
counterpart. Further, the mutagenic manipulations required to observe crypt 
heterogeneity may alter the normal behaviour of stem cells. 

As the above observations demonstrate, the intestinal crypt is a complex, 
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highly regulated system whose dynamics under normal and pathological con¬ 
ditions are difficult to interpret using verbal arguments alone. Mathematical 
modelling offers a powerful tool that can provide mechanistic insights that com¬ 
plement and reinforce knowledge acquired from molecular, histological and live- 
imaging studies. Models can be used to integrate our assumptions and measure¬ 
ments within a consistent theoretical framework, test competing hypotheses in 
silica, and generate new predictions that can then be validated experimentally. 
In the next section we provide a brief overview of mathematical approaches 
that have been used to describe crypt dynamics, before focussing in Section 0] 
on multiscale models of crypt organization and carcinogenesis. 


3 Modelling approaches for colorectal cancer 

Mathematical modelling has been used to investigate aspects of colorectal cancer 
for over half a century. The first such model, by Armitage and Doll, [3] was used 
to explain experimental data indicating a power law relationship between age 
and cancer incidence. This model highlights the multistage nature of colorectal 
cancer, by predicting that approximately six successive mutations in a single 
cell are required for carcinogenesis to occur. 

As the quality of cellular-resolution experimental data has improved over the 
last two decades, a variety of mathematical models have been used to study cell 
population dynamics in the colonic crypt. These include both deterministic [TUI 
[36] and stochastic [44l [56] models. As an illustrative example we highlight the 
work of Johnston et al.|5B], who analyse and compare age-structured and con¬ 
tinuous models of crypt population dynamics and propose two feedback mech¬ 
anisms that may regulate cell numbers to maintain crypt homeostasis. Here 
we focus on their continuous model, which considers timescales much greater 
than the typical crypt cell cycle time and populations large enough to warrant 
a continuous approach. Letting Ab(t), Ni{t) and N 2 {t) denote the size of the 

e set of coupled ordinary differential equations 

(as — «! — a2)No, (1) 

(/33-/3i-/32)iVi+a2iVo, (2) 

P 2 N 1 -JN 2 . (3) 

id as (respectively /3i, /32 and /Ss) denote the per- 
capita rates of stem (respectively transit amplifying) cell death, differentiation 
and proliferation, while 7 denotes the per-capita rate of removal of differentiated 
cells due to sloughing. If these rates are assumed to take constant values, then 
given an initial crypt population Nq{0) = uq, fVi(O) = hi, iV 2 ( 0 ) = h 2 , the 


stem, this model comprises 

dfVo 

dt 

dNi 

dt 

dN2 

dt 

Here the parameters ai , a 2 
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system of equations (Cl)-® has the straightforward solution 


Noit) = noe“*, 

(4) 

Afi(t)=Ae“*-k(ni-A)e^‘, 

(5) 

N2{t) = -k Ce^* + {n2 - B - C)e-'^\ 

(6) 


where a = 03 — oi — 02 , /3 = — /3i — /32 and 


Oi2no P 2 A 

n ^ ^ I 5 

a — p 7 + a 


and 


_ /32(hi - A) 

1 + P 


(7) 


This solution exhibits exponential growth or decay according to the values of a, 
P and 7 . Johnston et al.|36j note that such structural instability is unrealistic, 
since a crypt must be able to maintain homeostasis in the face of fluctuations in 
rates of proliferation and differentiation, for example due to stochasticity. Some 
form of feedback is required to ensure such robustness. The authors examine two 
possibilities in their model, feedback on proportions of cells differentiating and 
on rates of cell division, and deduce that the latter form offers a more feasible 
explanation for observed lag phases in crypt dynamics after mutations occur, 
and thus the existence of benign tumours prior to the onset of carcinogenesis. 
Such work highlights the utility of simple, analytically tractable population 
dynamics models in this field. 

While earlier models focus on temporal dynamics only, increasingly mod¬ 
ellers have sought also to understand the spatial organisation of the crypt. The 
classical approach to modelling spatial dynamics is to treat the tissue as a 
continuum. [24] As such models average over length scales that are much larger 
than the typical diameter of a cell, processes such as proliferation and adhesion, 
which may become disrupted in cancerous cells, can be considered only in an 
averaged sense. Such models do not allow the straightforward incorporation of 
subcellular processes, such as signalling pathways or protein-level descriptions of 
cell-cell adhesion. When studying colonic crypt organisation and carcinogenesis, 
a further challenge to the use of classical, continuum models is the relatively 
small numbers of cells in the tissue of interest; intestinal crypts contain around 
300 cells in mice and 2000 cells in humans. 

As a result of these limitations, over the last few years a number of authors 
have adopted a discrete approach to modelling crypt cell dynamics in health 
and disease. Several discrete approaches have been applied to model crypt cell 
population dynamics and combine descriptions of intracellular and cellular pro¬ 
cesses. These vary in complexity from fixed-lattice cellular automata and the 
cellular Potts model to off-lattice cell-centre and vertex models. Such multiscale 
models allow for the integration of disparate experimental data and the predic¬ 
tion of system-level behaviour, but are more complicated and less amenable to 
mathematical analysis than their single-scale counterparts. pS] 

Before proceeding to focus on multiscale modelling, we note the availability 
of recent review articles on other aspects of intestinal tissue renewal and col¬ 
orectal cancer. These include reviews focused on the biological insights gained 
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by compartment models of crypt dynamics, m the contrasting approaches and 
software implementations for spatial computational models. [T9] and the bur¬ 
geoning links between modelling and experiments in this field. [38] 


4 Multiscale models of crypt organization and 
carcinogenesis 

We next focus on a small number of representative individual cell-based and 
multiscale models of crypt dynamics, broadly in order of mathematical com¬ 
plexity. 

4.1 Stem cell dynamics in the crypt 

One of the first theoretical crypt studies to explicitly model the spatiotemporal 
dynamics of individual cells is by Loeffler and coworkers. [43] who use a stochastic 
cellular automaton approach. In cellular automata, physical space is discretized 
using a fixed lattice on which cells can occupy individual sites. The rules for 
cell movement are formulated in terms of cells moving between lattice sites. 
The state of each cell in the system is updated over discrete time steps using 
automaton rules and the state of the cell’s neighbourhood at the previous time 
step. [U Ho] It is possible to efficiently simulate the interactions of large numbers 
of cells in cellular automata; this advantage comes at the expense of an artificial 
spatial lattice anisotropy and hence discontinuous cell motion. 

In their model, Loeffler and colleagues [231 [SH] demonstrate that the spatial 
arrangement of proliferating columnar cells, as well as mucus-producing Goblet 
cells, can be explained by the concept of a cellular pedigree. This assumes that 
a self-maintaining stem cell generates transit amplifying cells, which undergo 
a well-defined series of cell divisions before ceasing proliferation. Under this 
hypothesis, regulatory control is required only in the stem cell region. Loeffler 
et al.|23] also investigate how the probability of a stem cell generating zero, one 
or two stem cells upon division affects the crypt structure. The authors initially 
assume these probabilities to take constant values, thus neglecting possible feed¬ 
back mechanisms, and find that the probability of asymmetric division should 
lie in the range (0.8, 0.95), assuming a 24 h cell cycle time. This prediction 
proved to be inconsistent with more recent data on the time scale over which 
mutated crypts become monoclonal, m highlighting how such models may be 
validated as new measurements become available. 

In further work, Loeffler and colleagues seek to modify their crypt model 
so as to recapitulate the observed timescales of monoclonal conversion. [S^ The 
authors incorporate a state dependence into the probabilities of stem cell divi¬ 
sion and stem cell cycle time, in which stem cell growth is inhibited if a crypt 
already contains numerous stem cells. Crypt fission is assumed to occur at 
some threshold number of stem cells per crypt. Their modified model is shown 
to be consistent with data on crypt fission and monoclonal conversion times. 
This research demonstrates that some form of autoregulatory feedback is likely 
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to modulate rates and modes of stem cell division in response to fluctuating 
environmental conditions, reflected for example by varying crypt cell numbers. 


4.2 Crypt homeostasis and monoclonal conversion 

The first off-lattice model of crypt dynamics was developed by Meineke et al. [48] 
We choose to focus on the mathematical details of this model as it forms the 
basis of several more recent crypt models. [Tiisniiiiiiiiiii] Unlike lattice-based 
cellular automata, off-lattice models allow continuous cell positions and enable 
the study of mechanical effects on cell populations in a straightforward manner. 
Off-lattice models offer a more physically realistic view of cell population dy¬ 
namics, and may be parameterized by biophysical and kinetic parameters that 
can be determined experimentally. 

Meineke et al.jlS] model the crypt as a cylinder, which is unrolled to give 
a flat rectangular surface with periodic boundary conditions on the left and 
right sides. There are two basic components to the authors’ description of cell 
mechanical interactions. The first step is to decide which cells are neighbours 
of each other, the second is to determine the forces transmitted by these neigh¬ 
bours. Each cell i is represented by a point S which may be interpreted as 
the centre of the nucleus. Meineke et al.|48] assume that cell movement within 
the crypt is driven by mitotic activity. New-born cells force nearby cells to move 
away, causing a pressure-driven passive movement primarily up the crypt axis. 
The authors model the local repelling and attracting intercellular forces by a 
network of springs connecting neighbouring cells. The total force Fi(t) acting 
on a cell i at time t is equal to the sum of all forces coming from the springs of 
all neighbouring cells j € Mi{t) adjacent to i at that time, 

Xy (t) (||x„(t)|| - sy(t)), (8) 

jeMiit) 


where fi is the spring constant, x^ (t) = x.i{t) — Xj (t) is the vector from i to j at 
time t, Si-ij {t) is the corresponding unit vector and Sij is the natural separation 
(spring length) between cells i and j. For simplicity, all cells are assumed to 
have identical mechanical properties. 

Adopting the simplifying assumption that inertial terms are small compared 
to dissipative terms, Meineke et al.|48] arrive at first-order dynamics with the 
evolution of each cell determined by 


dxt 

dt 






(9) 


where ry is the damping constant. Equation (jH]) is solved numerically by using 
a simple forward Euler discretization, so that the effective displacement within 
a small time interval At is given by 


Xi(t + At) = Xi(t) -I- —Fi(t)At. 


( 10 ) 
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As this numerical method is explicit, a sufficiently small time step must be 
chosen. 

After any cell movement, the cell connectivity network is updated using a 
Delaunay triangulation and cell boundaries are determined by a Voronoi tes¬ 
sellation, which generates polygonal cell packing. This approach is motivated 
by experimental evidence that columnar cell boundaries have a convex polygo¬ 
nal shape of about equal area, as seen in projection from the intestinal lumen. 
In the model, stem cells are restricted to lie at the bottom of the crypt, and 
asymmetric stem cell division is assumed. 

Meineke et al.[lH] implement cell division in their model as follows. Upon 
division, each daughter cell is assigned a random cell cycle time drawn from a 
suitable distribution. When a cell is ready to divide, a new cell is introduced 
a short distance e away from the parent cell, at a random angle. To account 
for cell growth, and to avoid an unrealistically large force between these two 
neighbouring cells, the rest length of the spring between them is chosen to 
increase linearly from e to its normal value i over the first hour after division, 


r e-b (^ - e)(t - id), t<td + l, 

t > id -b 1, 


( 11 ) 


where td denotes the time at which the parent cell divides. 

Considering cell proliferation, Meineke et al.[48] follow earlier work[4^ by 
employing the cellular pedigree concept. Their off-lattice model reproduces data 
from crypt labelling experiments as well as the cellular automaton model of Loef- 
fler et al. [131 [SI] without the need for the state-dependent cell cycle assumption. 
Whilst resolving a number of issues with previous modelling frameworks, the 
model by Meineke et al.|48] is unable to reproduce the broad, wavy ribbons of 
clonal populations observed experimentally. m This model limitation is a con¬ 
sequence of the assumption that stem cell positions are fixed, which enforces 
asymmetric division and thereby prevents propagation of a labelling mutation 
within the stem-cell niche. 

To address this and other issues, van Leeuwen et al. ES] extend the model by 
considering the role of Wnt signalling in the crypt. The inclusion of a Wnt gradi¬ 
ent along the crypt axis that defines the stem cell niche allowed the assumption 
that stem cell positions are fixed to be relaxed. The local Wnt concentration 
is assumed to determine intracellular gene expression levels affecting the cell 
cycle, with high levels of Wnt at the crypt base driving cell proliferation and 
defining the stem cell niche. A threshold of Wnt signalling is assumed to be 
necessary for mitosis to occur. This crypt model represents the first attempt 
to link phenomena at the subcellular, cellular, and tissue levels of organization. 
In the model, Wnt signalling is assumed to control cell cyle progression, as well 
as mechanical interactions through the incorporation of contact edge-dependent 
cell-cell adhesion and cell size-dependent cell-matrix adhesion. The model re¬ 
capitulates the observed spatial variation in nuclear /3-catenin levels along the 
crypt axis. Notably, the model exhibits a gradual change in cell cycle times 
along the crypt axis, which suggests that, in contrast to the simple classifi- 
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Figure 2: Evolution of a mutant in the crypt model by Mirams et al.[50], based 
on earlier work by Meineke et al.[35] and van Leeuwen et al.[75]. A cell with 
a neutral mutation is introduced near the base of a homeostatic crypt and 
its progeny are tracked over time. Simulation snapshots are shown going left to 
right then top to bottom. Non-mutant cells are coloured yellow or red according 
to whether they are proliferating or differentiated, while all mutant cells are 
coloured blue. Image adapted with permission from Mirams et al. [50] ©Elsevier. 


cation of ‘stem’ and ‘non-stem’, cells may in reality adopt a more continuous 
representation of differentiation. 

By relaxing the assumption that stem cell positions are fixed, the van Leeuwen 
et al.[76] model recapitulates the observed processes of clonal expansion and 
monoclonal conversion (Eig. [5]). This aspect of crypt behaviour was studied in 
more detail by Mirams et al. , who show theoretically that the probability of 
crypt domination by a mutant population is increased if Wnt signalling is raised 
near the base of the crypt. This emergent property incorporating wild type and 
mutant cells, different levels of Wnt signalling, proliferation, and adhesion of¬ 
fers a possible explanation for the experimental observation that basally located 
stem cells are preferential oncogenic targets. [7] 

In the model by van Leeuwen et al.[76], the assumption of a cylindrical 
geometry can lead to unrealistic cell dynamics, especially near the bottom of 
the crypt. In more recent work, Fletcher et al.[25] guarantee a small number of 
stem cells at the bottom of the crypt by utilising a more realistic crypt geometry. 
This is achieved by projecting the crypt surface, approximated by a surface of 
revolution 


j: = {{r,e,z)eR^ : z = f{r)}, (12) 

onto the plane 2 = 0. This allows the planar model to be used still, facilitating 
efficient simulation, but with a different metric, determined by the projection. 
In more detail, let F^ define the spring force between two cells i and j given in 
the summand in equation ([5]), and let this force act in three-dimensional space. 
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Figure 3: Schematic illustrating the force calculation in a plane projection model 
of the crypt developed by Fletcher et al.[25]. See main text for details. 


Further, let Ilj and Uj respectively denote the tangent plane and outward unit 
normal to S at the point Xj. Constraining cells to lie on E results in a reaction 
force Rj acting in the direction —fij, whose magnitude balances the component 
of the force normal to Hj, denoted Fb. The component of the force F^ 

parallel to Ilj, denoted Fj'^, can be found by taking cross products. Using 
Lagrange’s formula for triple vector products, and letting denote the unit 
vector in the z direction, Fletcher et al.[25] arrive at the following expression 
for the projection of Fy onto the plane z = 0: 

Fy = Fij — {Fij ■ nj)hj — {Fij ■ e^)ej, + (F^- • nj)(nj ■ (13) 

A diagram illustrating this calculation is shown in Fig. [S] 

The geometrical approximation adopted by Fletcher et al.[25] is valid as long 
as the points are sufhciently close, which is likely the case for packed cells, and 
introduces a small error which increases with height up the crypt. This is in 
contrast to the cylindrical crypt model, giES] in which the error in approxi¬ 
mating the shape of the crypt is greatest at the bottom, the region which is 
critical to the process of monoclonal conversion. Using this model, Fletcher et 
al.[25] investigate different hypothesised modes of stem cell division by com¬ 
paring the typical time scale for monoclonal conversion under each hypothesis. 
The authors find that the probability and timescale for this process to occur are 
sensitive to the precise geometry of the stem cell niche and to the Wnt stimulus 
threshold necessary for mitosis to occur. 
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4.3 Cell fate determination and lineage specification 

The previous examples highlight how multiscale modelling has yielded mech¬ 
anistic insight into crypt stem cell dynamics and monoclonal conversion, and 
how these processes may be influenced by intracellular processes relating to the 
Wnt pathway. However, cell proliferation and differentiation in the crypt is not 
dependent solely on the presence of a Wnt ligand gradient. Other signals, such 
as Notch, may act in combination with Wnt signalling in this context. HHEg 
Recent crypt models have sought to understand how the interplay of these sig¬ 
nals regulates cell proliferation. In particular, following an approach similar to 
that of van Leeuwen et al.[75], Buske et al.[T2] study the combined influence of 
Wnt and Notch signalling on cell proliferation and differentiation within the first 
three-dimensional cell-based model of the crypt. Notably, the authors also pro¬ 
pose the first explicit representation of the basement membrane as a provider 
of mechanical integrity for the crypt, describing this by a three-dimensional 
triangular mesh of elastic fibres. In contrast to the Delaunay description of 
cell neighbours employed by Meineke et al. [3S] , the authors use an ‘overlapping 
spheres’ approach. Here, the radius Ri of each cell is accounted for in addition 
to its position x^. By considering changes to cell adhesion, deformation and 
compression energy terms arising from cell movements and interactions, the 
authors arrive at a set of coupled ordinary differential equations of the form 



(14) 


(15) 


where the coefficients tibm, VvOi VC describe friction between a cell and the 
basement membrane, in response to volumetric changes and between two cells, 
respectively. Here Afj denotes the contact area between two overlapping spher¬ 
ical cells i, j, while and Gi denote generalized forces acting on cell i. It is 
interesting to compare this more detailed biophysical description of cell move¬ 
ment with the simpler spring-based approach discussed above; we return to the 
issue of model comparison in Section 15.31 

Buske et al.|12j compute cellular Notch activity by summing contributions 
from cell neighbours and determining cell fate from a ‘look-up table’ that spec¬ 
ifies distinct outcomes for particular Wnt and Notch levels. This model re¬ 
capitulates experimental observations on the spatial distributions of secretory 
Paneth and Goblet cells in wild-type crypts [?,[59], timescales for monoclonal 
conversion [4T] and the effect of Wnt/Notch mutations on crypt organization. p51 
usiinaiis] as illustrated in Figure 01 Other models that have sought to address 
the rules underlying cell fate determination and lineage specification include the 
work of Pin et al.[5n] and Wong et al. m These models largely restrict their 
focus to simple rule sets for cell regulation, rather than extending the more 
detailed ordinary differential equation models of Wnt signalling (42] to account 
for crosstalk between pathways. 
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Wild type Wnt+ Wnt- Notch+ Notch- 



Figure 4: Simulation results for the model of cell fate determination in the crypt 
by Buske et al.[T2], showing the spatial distribution of undifferentiated/stem 
(red), Paneth (green), Goblet (yellow) and enterocyte progenitor cells (blue). 
Alongside the wild type case, four gain or loss of function perturbations are 
shown, each of which agree qualitatively with experimental observations. From 
left to right, these are: wild type; constitutive activation of Wnt in all cells 
(Wnt-I-), which results in loss of Goblet and enterocyte progenitors [65]: reduced 
Wnt signalling (Wnt-), which results in loss of undifferentiated and Paneth 
cells[23]; constitutive activation of Notch (Notch-1-), which results in loss of se¬ 
cretory lineages [2^: and inhibition of Notch signalling (Notch-), which results in 
loss of undifferentiated and absorptive cells. m Image adapted with permission 
from Buske et al.[12] under a Greative Gommons License. 
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4.4 Crypt deformations 

We conclude this section by highlighting the more recent application of discrete 
and multiscale modelling to study later steps in colorectal carcinogenesis. Fol¬ 
lowing the accumulation of genetic mutations which disrupt the dynamics of a 
crypt, excess proliferation can cause a build-up of stress, which may eventually 
cause the crypt to buckle and deform. Drasdo and Loeffler[20] were the first to 
investigate this process at a cellular level by modelling a vertical cross-section 
through a crypt as a U-shaped chain of growing, deformable elastic spheres. The 
authors demonstrate that crypt buckling can occur when cell cycle durations 
are reduced or the rigidity of cells is reduced. More recent work has sought 
to address the dynamic role of the basement membrane in determining crypt 
geometry [ 2 TI122 ) and the biomechanical behaviour of in vitro crypt organoid 
cultures, which reproduce the coarse three-dimensional structure of crypts and 
allow for more precise experimental measurement, imaging and perturbation. m 


5 Mathematical challenges 

Having summarised recent efforts to apply multiscale modelling techniques to 
study colonic crypt organization and colorectal carcinogenesis, we next review 
some of the key mathematical challenges in this area. 

5.1 Inference 

An important limitation associated with multiscale models of tissue dynamics 
is the lack of established statistical and computational techniques with which 
they may be parameterised or validated against experimental data. It may be 
argued that to date, the majority of multiscale models of carcinogenesis have 
been used to obtain qualitative insights into the spread and fixation of muta¬ 
tions within epithelial tissues. However, some recent studies have attempted to 
place such models on a more quantitative footing. Further, a variety of inverse 
problem frameworks have been developed for the estimation of model param¬ 
eters. Key issues to be addressed in this area include how to perform model 
selection, [61] develop simplified ‘emulators’ to approximate such models for pa¬ 
rameter inference. |33j and efficiently compute statistical likelihoods or perform 
approximate Bayesian computation (ABC). [80] 

While these issues have been actively pursued in the field of mathemati¬ 
cal ecology, to date relatively few studies have considered them in the context 
of tissue homeostasis and carcinogenesis. Of particular relevance is the work 
of Sottriva and Tavare, m who demonstrate how ABC may be applied using 
methylation pattern data [81] to infer posterior distributions for parameters of 
stem cell dynamics in a cellular Potts model of a colonic crypt. This model 
shares many features with the models discussed in Section [4] 
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5.2 Computational considerations 

A further barrier to the wider use of multiscale models for the study of tissue 
homeostasis and carcinogenesis is the lack of standards or benchmarks. Despite 
their establishment in the cancer modelling literature, it is still the case that pre¬ 
viously developed models and methods are seldom re-used effectively, because 
they are typically not available as rigorously tested, open-source simulation soft¬ 
ware. It is therefore difficult to guarantee the reproducibility of computational 
results, which is of particular importance given the complexity of such models. 

This problem is beginning to be addressed, for example by the Chaste 
project, which provides an open-source C-|—I- implementation of several classes 
of individual cell-based and multiscale models of multicellular populations. m 
Chaste constitutes a multi-purpose software library that is developed to sup¬ 
port computational simulations for a wide range of biological problems. It is 
developed using an agile approach and comprises fully tested, industrial-grade 
software. m A major component of this software library is the code used to sim¬ 
ulate several of the multiscale crypt models discussed in Section H) Full details 
of all technical aspects of the implementation of such models, and the code, are 
freely available to the community (http://www.cs.ox.ac.uk/chaste) along 
with documentation and tutorials. However, a more systematic and compre¬ 
hensive description of specific multiscale models in the literature remains an 
ongoing challenge. 

5.3 Effect of modelling framework 

To date, there has been little comparative study of different cell-based modelling 
approaches. There are clear strengths and weaknesses associated with each ap¬ 
proach, and therefore some situations in which it is clear which approach is 
most valid; for example, cellular automata are not well suited to the detailed 
study of cell adhesion and motility. However, there are other cases where it 
is not clear which modelling approach is valid. This raises the question of to 
what extent are the results of model simulations artifacts of the chosen mod¬ 
elling approach or method of numerical solution? It can be difficult to make an 
accurate comparison between different models, in order to evaluate the impact 
of various constitutive assumptions on the behaviour of a model. This difficulty 
is particularly relevant in the context of multiscale models, which may couple 
the biomechanical behaviour of cells in tissues with other processes such as the 
secretion, transport and uptake of nutrients or signalling molecules. When com¬ 
paring different constitutive assumptions using the same overall modelling ap¬ 
proach, a computational framework is required that allows one to easily change 
the fine details of a model and its implementation. A small number of studies 
have begun to address these issues in the context of crypt dynamics and ep¬ 
ithelial tissue dynamics more generally. [271 El] Other work has also sought to 
characterise the mechanical behaviour of distinct classes of individual cell-based 
models under loading, unloading and shearing. [HI [58] Such studies could help 
refine the region of parameter space that must be searched when fitting such 
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models to data. 


5.4 Coarse-graining approaches 

As the above examples demonstrate, multiscale modelling has played an increas¬ 
ingly important role in the study of colonic crypt organization and colorectal 
carcinogenesis. A key strength of this approach is that it allows the simulation of 
observed cell-level behaviours, such as cell heterogeneity and the occurrence of 
mutated cells. However, the computational cost associated with analysing such 
models remains a significant challenge. This is exacerbated by the stochastic 
behaviour of cell processes such as proliferation and differentiation, as well as 
by the increasingly detailed mathematical descriptions of cell decision-making 
processes. For example, in order to build up statistical distributions for how the 
probability and timescale of mutant domination varied with properties such as 
the location of the initial mutation in a colonic crypt, Mirams et al.|50) needed 
to run a large number of realisations of their model, cumulatively equivalent to 
over 9000 years of simulated ‘crypt time’. 

One approach to tackling this problem is to develop coarse-grained models 
that capture the essential features of the original models but are amenable 
to efficient simulation methods or established mathematical techniques (e.g. 
asymptotic or bifurcation analysis). As the methodology used to derive coarse¬ 
grained equations is strongly dependent on the class of underlying cell-based 
model, here we highlight a number of cases that are particularly relevant to the 
type of crowded cellular environment found along the crypt axis. A common 
feature shared by the coarse-grained models considered below is that, in one 
spatial dimension, the equation describing the spatiotemporal evolution of the 
coarse-grained cell number density q{y,t) takes the general form 



(16) 


Here y represents position along a spatial axis and D{q) is a nonlinear diffusion 
coefficient within which the details from the respective cell-based models appear. 

We start by considering work by Simpson and co-workers, who have ac¬ 
counted for volume exclusion effects in a study of stochastic cell movement on 
a lattice. [ll|3Zl|67] By way of example, we highlight the case of stochastic ex¬ 
cluded motion on a one-dimensional lattice. In this case, the authors find that 
the governing coarse-grained equation takes the form (1161) . with the nonlinear 
diffusion coefficient given by 


D{q) = Do{l - (Tg(4 - 3g)), 


(17) 


where tr is an adhesion parameter and Dq represents the diffusivity of an isolated 
cell. 

In other work, Alber and coworkers [H[5] show that the coarse-grained density 
derived from an underlying lattice-based cellular Potts model is also described 
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by m, with the nonlinear diffusion coefficient given by 


where K is a parameter that defines the maximum allowed cell density on the 
lattice and Dq is again the diffusion coefficient representing the random motion 
an isolated cell. Interestingly, as the density tends to the limiting value K, 
the diffusion coefficient tending to infinity represents the case of maximal cell 
packing. 

Considering off-lattice cell-based models, Murray and co-workers [52l [55] 
show how different inter-cellular force laws (e.g. linear. Hertz, Lennard-Jones) 
can be described using a range of corresponding nonlinear diffusion equations. 
As before, in one spatial dimension the governing equation takes the form (I16II 
and in the case of a linear spring force between cells, the diffusion coefficient 
takes the form 


m = ( 19 ) 

riq^ 

where k and rj denote the spring constant and damping coefficient, respectively. 
Subsequent studies have performed a similar derivation starting from alternative 
stochastic off-lattice cell-based models. ^51 

The work of Murray and coworkers [S2| |SS| IM] has been applied directly to 
the intestinal crypt models proposed by Meineke et al.[48] and van Leeuwen et 
al.[76] as described in Section [421 We therefore explore their model in further 
detail. The authors focus on the simulation framework described in Section [4.21 
a crypt of length L is populated by cells of equilibrium diameter a and treated 
as an unwrapped cylinder. Exploiting the axisymmetric nature of the crypt, 
the authors consider a model in one spatial dimension and show that the cell 
number density along the crypt axis, g(j/,t), satisfies the nonlinear diffusion 
equation 


where D{q) = k/ijq^, Tc is the mean cell cycle period in the proliferative region 
of the crypt, ycrit is the height up the crypt below which cells can proliferate 
and H(.) is the Heaviside function. In Figure [5] we compare measurements of 
the average cell velocity profile along the crypt axis from the underlying cell- 
based model with predictions from the continuum model for different points 
in parameter space. Despite the simplifying assumptions made in the deriva¬ 
tion of equation (1^01) . it can be seen that the coarse-grained model captures 
qualitatively the behaviour of the multi scale model. 

Upon deriving equation (1201) . it is natural to apply standard techniques of 
analysis. For example, a critical non-dimensional parameter is identified by the 
authors; seeking a steady-state solution in which the rate of proliferation in the 
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Figure 5: Steady-state velocity plotted as a function of crypt axis for different 
values of the spring constant, k, and cell cycle period Tc- Image reproduced 
with permission from Murray et al. [55] . Markers denote measurements from 
simulations of cell-based model. Lines denote solutions from continuum models. 
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base of the crypt is matched by the rate of loss at the top, the authors identify 
that the condition 


k > 


T] \'n‘2ycrit(,‘^kj ycrit) 


( 21 ) 


must hold in order that a steady-state solution exists. Whilst this result was 
verified using simulations of the underlying cell-based model, identifying such a 
parameter relationship using numerical simulations alone would be a non-trivial 
computational task. It is notable from the above examples that the nonlinear 
diffusion equations that have arisen from a number of disparate coarse-graining 
approaches provide a unifying framework with which to classify and charac¬ 
terise behaviour in the underlying cell-based model. A major mathematical 
challenge in this endeavour is to extend such work to more detailed two- and 
three-dimensional tissue geometries. 


6 Concluding remarks 

This review has highlighted the role that multiscale mathematical modelling 
can play in gaining mechanistic insights into the spatiotemporal organization 
of intestinal crypts under normal and pathological conditions. We have high¬ 
lighted major mathematical and computational challenges associated with the 
multiscale modelling approach, concluding with a discussion of recent efforts 
to obtain coarse-grained descriptions of such models. Such attempts offer one 
possible way of leveraging the well-developed tools of mathematical analysis for 
continuous models to address key problems in multiscale modelling. 
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